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ABSTRACT 

Periodic  Peltier  heating  and  cooling  may  be  generated  at  a  solid/ 
liquid  interface  by  the  passage  of  an  alternating  electric  current.  Malt¬ 
ing  and  freezing  will  then  occur  at  rates  which  will  be  influenced  by  the 
kinetic  law  of  molecular  attachment  relating  solidification  velocity  to 
interface  temperature.  The  equation  describing  the  interface  attach¬ 
ment  kinetics  may  be  determined  by  measuring  the  response  of  a  solid/ 
liquid  interface  to  square -wave  Peltier  heating  and  cooling.  The  math¬ 
ematical  methods  required  for  extracting  the  kinetic  equation  by  analy¬ 
sis  of  the  waveform  of  interface  motion  show  that  this  technique  should 
give  distinguishably  different  waveforms  for  linear,  parabolic,  and  ex¬ 
ponential  kinetics.  There  are  advantages  of  the  method,  but  precautions 
must  be  taken  to  eliminate  possible  confusing  effects. 


PROBLEM  STATUS 

This  is  a  final  report  on  one  phase  of  the  problem;  work  on  other 
phases  is  continuing 
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MEASUREMENT  OF  ATOMIC  KINETICS  OF  SOLIDIFICATION 
USING  PELTIER  HEATING  AND  COOLING 


INTRODUCTION 

At  present,  little  conclusive  information  is  available  concerning  the  atomic  mecha¬ 
nisms  by  which  metal  crystals  grow  from  their  melts.  Observations  of  morphological 
features,  such  as  crystallographic  facets,  have  usually  been  of  a  qualitative  nature.  An 
experimental  measurement  of  the  relation  between  crystal  growth  velocity  and  the  super¬ 
cooling  at  the  solid/liquid  interface  would  be  extremely  informative,  because  the  mathe¬ 
matical  form  of  this  relation  has  been  predicted  theoretically  on  the  basis  of  several 
proposed  crystal  growth  mechanisms.  The  present  report  describes  a  quantitative 
method  which  can  be  used  to  Investigate  this  relation  in  metals  and  semiconductors. 

Three  important  kinetic  relationships  have  been  predicted  on  the  basis  of  three  dis¬ 
tinct  atomic  crystallization  processes.  For  continuous  growth  at  atomically  rough  inter¬ 
faces,  where  all  surface  sites  are  equivalent,  the  Interface  velocity  is  expected  to  be  a 
linear  function  of  the  interface  supercooling  (1)  and  may  be  written 

V  =  (1) 

where  6  represents  the  temperature  deviation  of  the  interface  from  the  equilibrium  tem¬ 
perature,  and  nx  is  the  linear  kinetic  coefficient.  When  atoms  are  added  to  the  crystal 
at  surface  steps  generated  by  dislocations,  the  predicted  kinetic  relation  (2)  has  the  form 

<•  --  ^e\e\  .  (2) 

If  the  crystal  grows  by  the  independent  nucleation  of  successive  layers  on  an  atomically 
smooth  (singular)  surface,  an  exponential  kinetic  relation  (3)  is  obtained; 


The  absolute -value  signs  in  Eqs.  (2)  and  (3)  have  been  introduced  to  make  the  >-'s  odd 
functions  of  the  0’s,  on  the  assumption  that  freezing  and  melting  kinetics  are  the  same. 
While  the  measurement  of  an  experimental  relation  between  v  and  0  would,  in  principle, 
enable  one  to  d^tinguish  between  the  mechanisms  leading  to  Eqs.  (1)  through  (3),  in 
practice  great  difficulty  is  encountered  in  the  measurement  of  ?.  In  metals  substantial 
growth  velocities  are  expected  even  at  small  values  of  0,  but  such  rapid  growth  almost 
invariably  results  in  complex,  nonplanar  interfaces  with  large  temperature  gradients 
associated  with  the  release  of  the  latent  heat  of  fusion.  In  such  situations  the  tempera¬ 
ture  of  the  interface  cannot  be  measured  with  sufficient  accuracy  in  any  known  direct 
fashion. 

A  few  years  ago  Kramer  and  Tiller  (4)  proposed  an  experiment  li  which  interface 
velocities  and  temperature  variations  could  be  calculated  from  measurements  of  the  at¬ 
tenuation  of  a  thermal  wave  propagated  through  a  solid/llquld  interface.  The  great  sen¬ 
sitivity  in  that  experiment  was  made  possible  by  the  measurement  of  periodic  variations 
of  temperature  about  the  equilibrium  temperature,  rather  than  of  absolute  instantaneous 
or  steady-state  temperatures.  While  the  Krarrer-Tiller  thermal-wave  experiment  has 
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been  used  to  investigate  the  solidification  of  tin  (5,6)  and  ice  (7),  the  correctness  of  the 
execution  and  interpretation  of  that  experiment  has  been  questioned  (8-10).  Some  of  the 
more  serious  experimental  problems  have  been  the  prevention  of  convection  in  the  liquid 
and  the  elimination  of  "short  circuit"  paths  by  which  the  thermal  waves  flow  along  the 
container  walls  instead  of  through  the  solid/liquid  interface.  Also,  the  mathematical 
treatment  of  situations  in  which  the  kinetic  relation  between  interface  supercooling  and 
crystallization  velocity  is  nonlinear  has  been  very  difficult  (11).  Finally,  published  data 
on  tin,  derived  from  that  experiment,  have  been  conflicting  and  inconclusive  (5,6,8-10). 

We  now  describe  a  new  experiment  wnich  has  similarities  to  the  method  of  Krainrr 
and  Tiller  but  which  eliminates  some  of  the  experimental  and  analytical  difficulties  of 
their  approach. 


DESCRIPTION  OF  EXPERIMENT 

Our  experiment  derives  kinetic  data  from  an  analysis  oi  the  response  of  a  solid/ 
liquid  interface  to  a  periodic  input  of  heat  directly  at  the  interface;  this  heat  input  is  in* 
troduced  using  Peltier  (thermoelectric)  heating  and  cooling  generated  by  an  alternating 
electric  current  passing  through  the  Interface.  The  resulting  interface  motion  is  meas¬ 
ured  and  the  corresponding  interface  temperature  variation  calculated,  almost  the  re¬ 
verse  of  the  procedure  described  by  Kramer  and  Tiller,  who  measured  temperature 
variations  and  then  calculated  velocities. 

To  understand  the  response  of  a  solid/liquid  interface  to  Peltier  heating,  first  con¬ 
sider  a  system  solidifying  with  infinite  kinetics,  so  p  in  Eq.  (1)  is  infinite,  and  S  is  a!  • 
ways  zero.  When  a  direct  current  of  density  /  amp/cm2  passes  through  this  interface, 
Feltler  heat  is  emitted  or  absorbed,  depending  on  the  direction  of  current  flow,  in  the 
amount  <jP  -  r;0  per  unit  interface  area.  Here  /'  is  the  Pe'tier  coefficient  between  solid 
and  liquid.  Since  i  ,ie  assumption  of  infinite  kinetics  demands  that  the  interface  remain 
at  the  equilibrium  temperature,  latent  heat  must  be  generated  to  cancel  the  Peltier  1-sat 
exactly.  Thus,  the  interface  must  move  with  a  velocity  «0  =  /'/0  ky,  where  k  is  the  latent 
neat  of  fusion  In  joules  per  gram,  and  y  is  the  mass  density.  If  a  square -wave  electric 
current,  constant  in  magnitude  but  alternating  in  direction,  passes  through  the  solid/ 
liquid  Interface,  the  alternate  (Peltier)  heating  and  cooling  will  result  in  a  "sawtooth" 
motion  of  the  interface,  as  alternate  melting  and  freezing  occur  with  velocity  ±  v0 .  Joule 
heating,  which  is  Independent  of  current  direction,  will  remain  essentially  constant  with 
time  for  a  symmetrical  square  wave  and,  thus,  will  not  affect  the  periodic  motion  of  the 
solid/liquid  interface.  A  small  periodic  variation  in  Joule  heating,  due  to  the  change  of 
electrical  resistance  as  the  interface  moves,  Is  shown  in  a  later  section  to  be  negligible 
compared  to  the  periodic  variation  in  Peltier  heating. 

When  interface  kinetics  are  not  infinitely  rapid,  the  periodic  motion  of  the  interface 
must  be  associated  with  periodic  changes  In  Interface  temperature.  These  changes  in 
interface  femperature  imply  that  a  temperature  wave  is  propagated  into  the  solid  and 
liquid  material  on  either  side  of  the  interface  and  that  the  latent  heat  does  not  exactly 
cancel  the  Peltier  heat.  Under  these  conditions,  the  "sawtooth"  motion  of  the  solid/liquid 
interface  becomes  distorted,  and  If  this  distorted  waveform  can  be  measured,  then  the 
periodic  variation  of  t’.e  solid/liquid  interface  temperature  can  be  calculated. 

For  the  experimental  measurement  of  the  response  of  a  solid/liquid  interface  to 
square-wave  Peltier  heating,  the  material  to  be  investigated  is  enclosed  in  a  cylindrical 
container  equipped  with  four  electrodes  (Fig.  1).  This  container  is  positioned  in  a  tem¬ 
perature  gradient  so  that  the  lower  part  of  the  material  is  solid  and  the  remainder  Is 
liquid,  with  the  solid/liquid  interface  lying  between  the  second  and  third  electrodes.  The 
first  and  fourth  electrodes  introduce  the  square-wave  electric  current  which  generates 
Peltier  heating  and  cooling,  thus  causing  the  periodic  motioii  of  the  solid/liquid  interface. 
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Fig.  1  -  Apparatus  for  determination  of  solidifi¬ 
cation  kinetics,  showing  solid/liquid  interface, 
current  electrodes  (1  and  4),  and  potential  elec¬ 
trodes  ( 2  and  3) 


INTERFACF- 


S 


Since  solid  and  liquid  metal  have  different  electrical  resistivities,  motions  of  the  inter¬ 
face  cause  changes  in  the  resistance  of  the  metal  column  between  electrodes  2  and  3;  the 
interface  motions  are  detected  by  monitoring  these  resistance  changes. 

Considerations  of  Joule  heating,  radial  heat  flow,  and  sensitivity  of  measurement, 
indicate  that  the  optimum  diameter  of  the  metal  specimen  is  about  1  mm,  and  the  corre¬ 
sponding  amplitude  of  the  square-wave  current  is  in  the  range  of  2  to  5  amp.  Also,  it  is 
estimated  that  the  appropriate  fundamental  frequency  for  the  square-wave  current  induc¬ 
ing  the  interface  motion  should  be  roughly  10  Hz.  For  the  measurement  of  resistance 
variations  due  to  interface  motions,  a  high-frequency  (100  kHz)  sine-wave  current  ap¬ 
proximately  1  amp  in  amplitude  is  superimposed  on  the  square  wave.  Because  the  re¬ 
sistance  variations  due  to  interface  motion  are  small,  it  is  essential  that  a  bridge  method 
be  used  for  this  measurement.  Complete  details  of  the  electrical  circuitry  invoived  will 
be  presented  in  a  later  experimental  report. 

Experiments  of  the  type  under  consideration  are  applicable  only  to  substances  with 
sufficient  electrical  conductivity  to  allow  introduction  of  the  current  causing  the  Peltier 
heating  and  coolinf.  In  addition,  the  solid  and  liquid  must  differ  in  electrical  conductiv¬ 
ity  by  some  amount  sufficient  to  permit  accurate  resistometric  determinations  of  inter¬ 
face  motions.  Finally,  the  Peltier  c  ^efficient  at  the  solid/liquid  interface  must  be  large 
enough  to  produce  adequate  rates  of  heating  and  cooling.  Bismuth  and  germanium  show 
promising  combinations  of  these  properties,  and  several  other  elements  with  relatively 
low  melting  temperatures,  such  as  mercury,  gallium,  and  tin,  also  appear  to  be  suitable. 


ANALYSIS  OF  HEAT  FLOW 

We  will  now  show  how  interface  temperature  variations  can  be  calculated  from  meas¬ 
ured  waveforms  of  interface  motion.  We  shall  start  the  heat-flow  analysis  with  a  consid¬ 
eration  of  the  situation  in  which  an  interface  between  two  semi -infinite  substances  lies 
along  the  plane  x  o,  with  substance  1  occupying  the  space  j  >  o  and  substance  2  occupy¬ 
ing  the  space  x  <■  o.  One -dimensional  heat  flow  will  be  assumed  throughout  the  analysis. 
Let  heat  from  some  sinusoidal  source  with  angular  frequency  n  .*0  be  emitted  from  this 
interface  at  a  rate  of  yn  </n  sin  u  0/)  per  unit  area.  The  temperature  variations  re¬ 
sulting  in  each  substance  from  this  heating  will  be  given  by  solutions  of  the  heat-flow 
equation 


( 
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3^ 2 


with  the  boundary  condition 


o  . 


3-  >o 

» 

3-  <  0 


30 

i'/J  si»(nn  •  -*i  — 


+  *=  -r 


(4) 


(5) 


expressing  the  fact  that  no  heat  can  accumulale  at  the  interface.  The  upper  case  0  is 
used  to  express  temperature  as  a  function  of  x,  with  the  lower  case  8  always  indicating 
temperature  at  the  interfa  e  (« :  =  o)  only.  Here  the  a’s  and  k 's  represent  the  thermal 
diffusi'dties  and  conductivities,  respectively,  of  the  two  phases.  The  solutions  of  these 
equations  '’re  readily  obtained,  for  example,  by  an  extension  of  the  solution  given  by 
Jakob  (12)  and  may  be  written 


in  n - i 

y(  ncO 


_A _  -V(  n6j„  2  a  ] )  ; 


1 1n  I  sin  [n^t  -  naio  2a j)  x  -  v  '4] 


(6a) 


and 


A  +  VW„  -  .  , 


\TnaT) 


qn  \  sin  [iu0i  +  v(  «^0  -'J2)  *  ~  ”/*]  .  (6b) 


where  A  =  (v/v*7  +  *2  ’'oj)  1  •  Figure  2  illustrates  the  form  of  these  functions  at 

several  stages  during  the  cycle  for  a  wave  of  frequency  100  Hz  diffusing  from  a  boundary 
between  solid  and  liquid  bismuth.  Tt  is  seen  that  the  damping  is  extremely  rapid,  with 
the  oscillations  propagating  only  about  one  wavelength  into  either  medium  before  becom¬ 
ing  negligibly  small  compared  to  the  interface  temperature  variations. 


Fig.  2 -Thermal  wave,  frequency  100  Hz,  prop¬ 
agating  from  interface  between  solid  and  liquid 
bismuth.  The  wave  is  shown  at  three  stages: 

t2  =  t,  +  and  1 3  =  +  2w  3a>.  Dashed 

lines  show  envelopes  of  the  damped  sine  waves. 


The  expression  for  the  interface  temperature  variation  0n  can  be  obtained  by  set¬ 
ting  x  o  in  Eqs.  (6).  Hence, 


A  |  | 


n  xHTTl 


i7ni' 


vUV 


4) 


- 1  it/ 4 


(7) 
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where  we  have  used  a  complex  notation  for  the  sine  function  and  written  qn  for  the  sinus¬ 
oidally  varying  l?n  U‘  "a°  .  Wv  are  the  complex  notation  is  used,  it  is  to  be  understood 
that  the  actual  values  of  the  functions  represented  are  given  by  the  imaginary  part  of  the 
complex  functions. 

The  total  periodic  heat  emission  from  the  solid/liquid  interface  is,  in  fact,  the  alge¬ 
braic  sum  of  the  Peltier  and  latent  heat  emissions  and  may  be  expressed  as  the  sum  of 
sinusoidal  components,  so  we  may  write 


?(  t ) 


\<in\ 


i  (  n  a'0  t  ♦  <?n) 


qpUi  +g(U)  -  +\yv(t) 


(8) 


Since  qPU)  is  known  from  the  square-wave  electric  current  and  ?f(n  may  be  de¬ 
termined  by  measuring  »(<),  we  can  find  the  sets  of  |?J  and  <pn  which  describe  the  tota» 
interface  heat  emission.  The  corresponding  total  interface  temperature  variations  are 
then  calculated  from  Eq.  (7)  as 


«<•) 


I'/J 


i  (  nu0l  4) 

\n~ 


(9) 


The  calculated  0(t)  is  finally  plotted  against  the  measured  d(;)  to  reveal  the  experimen¬ 
tal  kinetic  relation. 


It  may  be  pointed  out  here  tl  at  this  experiment  uses  an  analysis  of  waveform, 
whereas  the  Kramer-Tiller  experiment  employs  measurements  of  the  amplitude  of  the 
transmitted  thermal  waves.  Those  authors  point  out  that  a  waveform  analysis  could  in 
principle  be  used  to  derive  the  kinetic  relation  from  their  experiment,  but  that  the  possi¬ 
bility  for  error  in  such  ar.  analysis  is  large.  The  present  experiment  would  measure 
waveforms  directly  at  the  interface,  rather  than  measuring  them  after  their  content  of 
higher  harmonics  has  been  reduced  by  conducting  through  some  length  ''f  material. 

Thus,  waveform  measurements  are  expected  to  be  more  accurate  n  the  present  experi¬ 
ment,  and  the  analysis  in  terms  of  waveform  more  straightforward. 


EXPECTED  WAVEFORMS 

It  is  desirable  to  have  quantitative  a  priori  knowledge  of  the  waveforms  of  interface 
motion  which  would  result  from  behavior  described  by  the  kinetic  relations  (1)  through 
(3).  These  waveforms  can  be  generated  by  either  a  Fourier-analysis  method  (for  linear 
kinetics)  or  an  integral -equation  method  (for  any  kinetic  law). 

To  solve  the  case  of  linear  kinetics,  we  first  note  that  each  frequency  component  of 
the  total  interface  temperature  variation  is  the  sum  of  Peltier  and  latent  heat  contribu¬ 
tions 


e 

n 


6Pn  +  °( 


(10) 


When  the  kinetics  are  linear,  we  can  express  the  latent  heat  emission  by  the  interface  as 
'Jin  -  ‘  K'>‘v „  '  ~^yp\e „•  11  we  now  consider  a  sinusoidal  interface  motion  due  to  a  sinusoi¬ 
dal  input  of  Peltier  heat  at  the  interface  with  frequency  n^0 ,  we  have,  from  Eq.  (7), 
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where  we  have  used  B  -  Substituting  Eq.  (11)  into  Eq.  (10),  we  can  solve  for 

t‘n  as  a  function  of  6Pn’  obtaining 

e  r  --  —  (i  .  JLe-'*  4  f  6;  e-2*" ■* 

"  1  +  ( /  4  n2)  \  v'T  n 

T’us  expression  relc*es  the  total  interface  temperature  variation  to  the  temperature 
variation  due  to  Peltier  heat  alone;  the  interface  velocity  is  then  given  by  Eq.  (1).  This 
velocity  may  be  related  to  the  original  sinus  lidal  input  of  Peltier  heat  by  expressing  0Pn 
in  Eq.  (12)  in  terms  of  qPr  by  the  use  of  Eq.  (7). 

For  the  square -wave  current  input  of  ampbtude  ;o  amp/cm2  and  angular  fundamen¬ 
tal  frequency  t*>0,  we  can  express  the  Peltier  heating  as  <lp  -  >-qPn,  with 


.  J*”  <  j  . 


(12) 


IPu  = 


e‘"^? 


n  odd. 


(13) 


The  interface  temperature  variation  can  be  written  e  =  16n,  with  the  en  values  given  by 
Eq.  (12). 


Equation  (7)  may  be  used  to  express  the  oPn  values  of  Eq.  (12)  in  terms  of  the  qPn 
valuer  of  Eq.  (13).  The  resulting  expression  for  P  may  be  normalized  by  multiplying  by 
the  kinetic  coefficient  to  give  the  velocity,  ar.d  dividing  by  the  velocity  for  infinite  ki¬ 
netics  t0.  The  normalized  velocity  Way)  may  be  written 


*0i ) 


\  n  ■ 

-  sin 

n2  +  B‘  L 


IT/  4) 


-  sin  (  -  Zv/\) 

\7> 


§1 

n 


mi  ( m<>Qt 


B3 

3v,'4)  -  —=■  sin  (n^,J  -  n) 


(14) 


Computer-generated  curves  of  this  function  are  shown  in  Fig.  3  for  various  values  of  b; 
for  B  <  io,  values  of  r.  up  to  199  were  used,  while  for  B=  lo ,  values  up  to  999  were  used. 
Also  shown  are  the  computer-generated  time  integrals  of  this  function,  giving  the  dis¬ 
placement  oj  the  interface  as  a  fraction  of  the  total  amplitude  the  oscillations  would  have 
with  infinite  kinetics. 

The  analysis  given  above  is  not  readily  extend'd  to  cases  of  nonlinear  kinetics. 
However,  an  alternative  method  employing  an  integral  equation  may  be  used,  if  we  make 
the  additional  assumption  that  the  interface  velocity  comes  close  to  the  steady-state  ve¬ 
locity  before  each  periodic  reversal  of  current  direction. 


In  this  method  we  consider  first  an  interface  between  two  semi-infinite  media  initially 
at  uniform  temperature  e  =  eo .  Suppose  that  at  time  t  -  t  '  this  interface  emits  a  small 
pulse  of  heat  of  magnitude  i>? .  The  temperature  distribution  in  the  two  media  at  times 
t  >  t  ’,  due  to  the  pulse  of  heat  at  time  t  -  t  ',  will  then  be  given  by 


a 

i 


SB . 


-  2  t  a.  (  l  - 1  ' ) 

\>(t  -  (T 


f  1.  X  >  0 

1  2.  I  <  0  ‘ 


(15) 


where  l  has  the  same  value  as  In  Eqs.  (6). 
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Fig.  3  -  Interface  displacement  and  velocity 
waveforms,  generated  from  Eq.  (14)  for  in¬ 
terfaces  with  linear  kinetics.  The  curves  have 
been  normalized  by  dividing  by  the  displace¬ 
ment  and  velocity  of  an  interface  with  infinite 
kinetics. 


The  interface  temperature  is  given  by  setting  x  =  o  in  Eq.  (15): 

e  =  et  +  &0=e  +i >? —(< -<  V'2.  (16) 

If,  instead  of  emitting  a  single  pulse  of  heat,  the  interface  emitted  heat  continuously  at  a 
rate  of  ?(«'),  starting  at  t  '  =  n ,  we  can  write  s?  -?(/')  dt  ' ,  and  the  total  temperature 
at  the  interface  will  be  given  by 


P  r  Q 

0 


) 


<tl  '  . 


(17) 


If  now  an  interface  is  initially  moving  at  its  steady -state  velocity  ur.-jr  the  influence 
of  a  direct  current,  ‘he  Peltier  heat  is  canceled  by  the  latent  heat  and  no  net  heat  is 
emitted  by  the  interface;  the  temperature  in  both  media  is  then  uniform  at  .  If  the 
current  is  reversed  at  time  t  =o,  the  function  <?(  t  ')  in  Eq.  (17)  will  be  composed  of  two 
terms,  due  to  Peltier  and  latent  heat,  now  no  longer  canceling: 

git')  =  Pja  -  Kyvieu  ')  j  .  (18) 


Substituting  into  Eq.  (17)  and  performing  the  integration  on  the  first  term,  we  obtain 


fA()  -  f50  +  ~r\2Pj0  yr  -  Ky 

V77 


f'  viou  ')  ] 

J,  -f) 


dt‘\  • 


(19) 


Substituting  the  kinetic  relations  (1)  through  (3)  into  this  expression  yields 


( 
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OU)  --  <?0 


and 


«*»  =  fc'o 

wnere  Eq.  (20a)  has  been  normalized  the  same  way  Eq.  (14)  was;  similar  normalizations 
are  not  useful  for  the  nonlinear  kinetic  laws.  The  values  of  eo  in  Eqs.  (20b)  and  (20c) 
are  determined  from  the  kinetic  laws  (2)  and  (3),  assuming  the  interface  to  be  mov^g 
initially  at  the  steady -state  velocity. 

Equation  (20a)  may  be  solved  by  applying  a  Laplace  transformation,  yielding 

W-V>  =  7-  =  2fSJV  erfc  (Bx'V)  -  1  .  (21) 

0 

but  the  nonlinear  nature  of  Eqs.  (20b)  and  (20c)  precludes  their  solution  by  this  method. 
Instead,  numerical  solutions  to  Eqs.  (20)  may  be  obtained  by  replacing  the  integrals  with 
series  approximations,  and  computing  the  interface  temperature  after  successive  time 
intervals  A  / ; 


+  ~P>0 


-  _[_(■'  0(t  ')  |6Hf  ')  \dt  ' 

V  '-V-U.  ^TT‘. 


(20b) 


+  4'''° 

V71 


2  vT 


-  C 


Jt  *  =0 


e( 


ut 


I eu  ')  |  y/u  -t‘) 


(20c) 


\  (nAo>0O  r  1 - y/(&UJQt ) 


*  ■  \/(  - 


V(  »  “  m) 


(22a) 


£„UAO  =  e0  1 


2  V" 


-  1 


0  (to  V(  n  -  m) 


(22b) 


and 


9„(nAo  =  6>0  ♦  -r:Pi0  \ (TST) 

,r»7 


m=o  -  m) 


(z2c) 


Figure  4  compares,  on  an  expanded  scale,  the  function  Way)  as  calculated  from  (a) 
the  Fourier  analysis,  Eq.  (14),  (b)  the  analytical  solution,  Eq.  (21),  of  the  integral  equa¬ 
tion,  Eq.  (20a),  and  (c)  the  series  approximation,  Eq.  (22a).  The  discrepancy  between 
curves  a  and  b  is  due  to  the  approximation  in  the  integral -equation  method  that  the  in  - 
terface  actually  attains  its  steady-state  velocity  and  temperature  before  the  current  re¬ 
versal;  the  discrepancy  between  curves  b  and  c  represents  the  error  due  to  approxi¬ 
mating  the  integral  by  the  series.  Since  these  discrepancies  are  small,  it  is  concluded 
that  the  functions  t-„(  vV )  computed  from  Eqs.  (22b)  and  (22c)  will  be  valid  representa¬ 
tions  of  the  temperature  changes  of  solid/liquid  interfaces  with  parabolic  and  er.por.ential 
kinetics  when  subjected  to  a  single  current  reversal.  The  response  of  these  interfaces 
to  a  square-wave  current  may  be  approximated  by  plo*tinu  several  reversals  in  succes¬ 
sion.  The  velocities  of  the  interfaces  may  be  computed  from  the  temperatures  by  use  of 
the  kinetic  relations,  Eqs.  (1)  through  (3);  r(t )  ,  the  position  of  the  interfaces,  may  be 
computed  as 
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Fig.  4  -  Velocity  waveform  for  interface  displaying  linear  kinetics 
with  B  -  10,  calculateJ  by  three  different  methods:  (a)  Fourier 
analysis,  Eq.  (14);  (b)  analytical  solution,  Eq.  (21),  to  integral 
equation;  and  (c)  numerical  computation  of  integral  equation 


For  linear  kinetics  we  obtain 


«••(/■)  i'. 


■  v0  '  p  e B  ,  0<  er  fc  (li  v'V  ’  -  >]  f  vT^H 


while  for  nonlinear  kinetics  the  integration  is  performed  numerically.  The  solution  for 
parabolic  kinetics  shown  in  Fig.  5  was  computed  using  At  =  10  1  sec,  APj  0  o.  l  /n,  and 
6  =  0.01  C°,  corresponding  to  m2  ~  3  cm/sec-JC2,  while  Fig.  6  shows  the  result  of  a 
computer  solution  to  Eq.  (22c)  for  a  steeply  varying  exponential  kinetic  law. 

The  computed  curves  indicate  that  when  the  solidification  kinetics  are  sufficiently 
sluggish  to  distort  the  waveform  of  interface  motion  from  a  simple  sawtooth,  the  differ¬ 
ent  kinetic  laws  (Eqs.  (1)  through  (3))  do,  in  fact,  produce  significantly  different  wave¬ 
forms.  Linear  kinetics  are  most  rapid  and,  therefore,  the  most  difficult  to  detect;  while 
the  values  of  m,  =  0.01  cm/sec-°C  reported  for  tin  by  Rigncy  and  Blakely  (6)  would  lead 
to  readily  observable  kinetic  effects  in  this  experiment,  these  values  appear  to  be  too 
low  to  be  valid  estim  ss  (7).  If  m,  is  greater  than  about  1  cm/sec-°C,  then  it  will  be 
extremely  difficult  to  measure  the  kinetics  by  means  of  this  experiment.  Parabolic  or 
exponential  kinetics  should  lead  to  more  easily  observable  interface  kinetic  effects. 
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Fig  5  -  Inte’  face  -displacement  waveform 
for  an  interface  displaying  parabolic  ki¬ 
netics  with  h-2  =  3  cm/sec-°C2,  subjected 
to  a  40-Hz  square  wave 


Fig.  6  -  Interface -displacement  waveform 
for  an  interface  displaying  steeplv  varying 
exponential  kinetics,  subjected  to  a  0.6-Hz 
square  wave.  For  this  interface,  ^  = 

1.4  x  10 19  cm  /sec  and  /it  =  50  C°  in  Eq.  (3). 


DISCUSSION 

Advantages  of  Peltier  Heating  Method 

The  use  of  the  Peltier  effect  to  introduce  heat  directly  at  the  solid/liquid  interface 
has  several  advantages  over  methods  in  which  thermal  waves  must  be  propagated  toward 
the  interface  through  the  specimen. 

The  maximum  amplitude  of  the  temperature  variation  occurs  directly  at  the  inter¬ 
face  and  is  equal  to  the  kinetic  supercooling.  Temperature  variations  diminish  rapidly 
as  one  moves  away  from  the  interface,  Eq.  (6),  so  that  frequency -associated  convection 
in  the  liquid  is  reduced  to  a  minimum.  The  "short  circuiting"  effect  of  periodic  heat  flow 
along  the  container  walls,  which  is  a  significant  problem  in  the  Kramer-Tiller  method, 
does  not  arise  here  since  no  transport  of  heat  to  a  temperature  sensor  is  required.  The 
waveform  and  amplitude  of  the  heat  input  at  the  interface  ?.r a  readily  controlled  by  the 
Peltier  effect;  the  use  of  a  square-wave  electric  current  is  a  convenient  way  to  keep  the 
Joule  heating  constani,  and  the  square  wave  is  a  useful  form  for  the  mathematical  analy¬ 
sis  of  the  response  of  systems  with  nonlinear  kinetics.  Finally,  direct  measurement  of 
the  waveform  of  interface  motion  by  an  electronic  method  here  makes  possible  the  de 
termination  in  a  single  experiment  of  the  f>'  tion  relating  interface  velocity  to  interface 
temperature. 


Justification  of  Approximations 

Several  effects  which  might  alter  the  analysis  of  heat  flow  as  presented  in  the  sec¬ 
tions  "Analysis  of  Heat  Flow"  and  "Expected  Waveforms"  have  been  neglected,  bu  argu  ¬ 
ments  are  presented  here  that  the  influence  of  these  effects  is  insignificant. 
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The  flow  of  part  of  the  periodic  heat  emission  from  the  Interface  out  ;  trough  the 
capillary  walls  would  invalidate  the  assumption  of  one-dimensional  heat  flow  and  thus 
affect  the  motion  of  the  interface;  this  effect  is  smaller  the  larger  the  diameter  of  the 
metal  specimen  and  the  higher  the  fundamental  frequency  of  the  square  wave.  However, 
increasing  the  diameter  of  the  specimen  would  require  large  increases  In  both  the 
square-wave  amplitude  and  the  higher  frequency  sine-wave  amplitude  If  the  sensitivity  of 
the  experiment  were  to  be  maintained;  moreover,  Joule  heating  would  then  produce 
troublesome  radial  temperature  gradients.  The  magnitude  of  the  electric  current  can 
be  kept  within  reasonable  limits  only  if  the  specimen  diameter  is  not  more  than  a  few 
millimeters.  The  range  of  useful  square-wave  (fundamentr.l)  frequencies  is  limited  by 
the  decrease  in  amplitude  of  interface  motions  at  high  freqjencles  and  the  approach  of 
the  Interface  motion  to  a  simple  sawtooth  at  low  frequencies.  Estimates  of  reasonable 
values  for  kinetic  coefficients  suggest  that  useful  kinetic  data  might  be  obtained  at 
square-wave  frequencies  of  roughly  10  to  20  Hz.  At  such  frequencies  the  decay  distance 
of  the  therm?'  waves  propagated  into  the  solid  and  liquid  metal  Is  short  [Eq.  (6)],  and  the 
area  over  which  the  container  walls  are  subjected  to  significant  changes  In  temperature 
is  small.  Since  the  thermal  conductivity  of  the  container  walls  Is  less  than  that  of  the 
metal  by  an  order  of  magnitude  or  more,  the  frequency -associated  heat  conduction 
through  the  container  walls  is  probably  unimportant,  even  In  a  1-mm  capillary. 

The  frequency-associated  variations  in  Joule  heating,  due  to  the  variation  in  speci¬ 
men  resistance  with  Interface  motion,  are  much  smaller  than  the  amplitude  of  the  Peltier 
heating  square  wave  Induced  at  the  Interface.  The  maximum  amplitude  of  the  Interface 
motion  (assuming  infinite  kinetics)  will  be  A'0  =  ’t»#/u>0  ”  woKy-  T*ie  maximum  vari¬ 
ation  In  Joule  heating  will,  therefore,  be  =  /’<  pt  -  ps)  A'#  ^r2  -  l2(pt  -  pa)  Pj0  r1, 
where  pt  and  p3  are  the  resistivities  of  the  liquid  and  the  solid,  and  l  Is  the  current. 
Since  the  Peltier  heat  peak-to-peak  amplitude  for  a  square  wave  is  A Qr  -2 Pi  ,  we  have 
AQj  &Qp  -  ( Pi  -  p3)  1 2  2kya >0r2  *  io  6  I-  for  a  10-Hz  square  wave  in  a  1-mm  capillary 
filled  with  bismuth.  Thus,  the  periodic  variation  In  Joule  heating  is  negligible  for  cur¬ 
rents  of  a  lew  amperes. 

Both  the  Fourier  method  and  the  integral-equation  method  treated  the  heat  f’ow  from 
a  stationary  Interface;  the  heat  is  considered  to  emanate  from  a  plane  fixed  in  space,  and 
the  change  of  position  of  this  plane  is  neglected,  although  the  latent  heat  resulting  from 
this  slight  change  of  position  Is  an  essential  part  of  the  analysis.  It  Is  believed,  however, 
that  the  error  Introduced  by  this  approximation  is  extremely  small,  as  the  amplitude  of 
the  interface  position  oscillations  Is  much  smaller  than  the  wavelength  of  the  thermal 
waves  emitted  by  the  interface.  From  Eq.  (6a)  we  can  see  that  the  wavelength  ol  thermal 
waves  of  angular  frequency  y  is  I .  =  2-r  %  2*  .  The  ratio  of  the  maximum  amplitude  of 

the  interface  motion,  given  in  the  preceding  paragraph,  to  the  thermal  wavelength  Is  thus 
A0  L  --  PjQ  /2>,  >■  \  2  aa  .  Substituting  typical  values  for  bismuth  of  P  =  0.02  v,  jg  =  100  amp/ 
cm2,  a  =  50  joules/g,  y  =  10  g/cm3,  and  *  =  0.05  cm\/sec,  we  have  A0  L  -  o.oos  v%- 
For  waves  with  fundamental  periods  of  a  few  seconds  or  less,  the  effect  of  interface  mo¬ 
tion  should  thus  be  very  small. 


Interpretation  of  Experimental  Results 

Thermal-wave  experiments  are  capable  of  determining  only  average  interface  veloc¬ 
ity  as  a  function  of  average  interface  supercooling.  O'Hara,  Tarshis,  and  Tiller  (13) 
have  suggested  that  if  some  crystallographic  facet  covers  a  significant  fraction,  but  not 
all,  of  the  interface,  then  a  combination  of  kinetic  processes  will  influence  the  experi¬ 
ment.  The  derived  kinetic  law  will  then  not  be  indicative  of  ?ny  single  kinetic  process. 

C  .  ra,  t  al.,  also  suggest  that  an  expe  rimental  finding  of  linear  kinetics  may  be  insuf¬ 
ficient  to  di  tinguish  between  a  truly  continuous  growth  process  and  a  process  involving 
growth  at  the  edge  of  layers  intersecting  the  solid/llquld  interface. 
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It  is  possible  that  a  thermal-wave  experiment  will  not  indicate  kinetics  described  by 
Sq.  (2)  when  growth  takes  place  by  a  screw-dislocation  mechanism.  Equation  (2)  was  de¬ 
rived  on  the  assumption  that  the  growth  layers  emanating  from  a  screw  dislocation  had 
assumed  their  steady-state  configuration,  so  that  the  spiral  growth  steps  rotated  on  the 
surface  without  changing  shape  (14).  If  the  direction  of  interface  motion  were  suddenly 
reversed,  the  entire  spiral  would  unwind  and  then  coil  up  in  the  opposite  sense.  During 
this  process  of  reversal,  however,  the  kinetics  would  not  be  described  by  Eq.  (2).  Only 
if  the  time  required  for  the  spiral  to  approach  its  steady-state  configuration  were  small 
compared  to  the  time  in  which  significant  changes  in  interface  temperature  occur  would 
the  system  behave  as  described  by  £q.  (20b). 

Let  us  consider  how  some  of  these  uncertainties  in  the  interpretation  of  experimen¬ 
tal  results  can  be  resolved.  Of  primary  importance  is  the  measurement  of  solidification 
kinetics  as  a  function  of  crystallographic  orientation.  Solidification  kinetics  as  described 
by  Eqs.  (2)  and  (3)  can  take  place  only  if  the  interface  lies  on  or  very  close  to  a  low -index 
crystallographic  plane  which  can  grow  only  by  the  generation  of  new  layers  by  screw  dis¬ 
locations  or  layer  nucleation.  Planes  which  are  far  from  low-index  planes  must  always 
contain  an  abundance  of  layer  edges  and,  therefore,  would  always  be  expected  to  show 
apparent  linear  kinetics  associated  with  layer  spreading.  So  long  as  these  layer  edges 
are  not  so  close  together  that  they  interact  with  one  another,  then  the  apparent  linear 
kinetic  coefficient  should  be  directly  proportional  to  the  density  of  layer  edges,  which  is 
readily  related  to  the  crystallographic  orientation  of  the  interface.  Therefore,  in  a  sub¬ 
stance  such  as  bismuth,  which  appears  both  macroscopically  and  microscopically  (15)  to 
grow  by  a  layer -spreading  mechanism,  it  is  necessary  to  study  a  variety  of  orientations 
near  those  which  show  facets.  For  example,  a  screw-dlslocatlon  mechanism  would  be 
Identified  if  we  found  parabolic  kinetics  at  and  very  near  the  low -index  plane,  and  (ap¬ 
parent)  linear  kinetics  with  orientation -dependent  coefficient  as  we  moved  away  from  this 
plane. 

Any  waveform  distortion  resulting  from  transient  coiling  effects  of  screw  disloca¬ 
tions  should  depend  on  the  amplitude  and  frequency  of  the  square  wave  driving  the  inter¬ 
face  motion;  we  can  conce  de  that  these  effects  are  not  important  If  we  obtain  consistent 
results  at  different  amplitudes  and  frequencies. 


CONCLUSIONS 

The  atomic  kinetics  of  solidification  of  some  substances  can  be  investigated  by  an 
analysis  of  the  response  of  these  systems  to  a  square  wave  of  Peltier  heat  input  at  the 
solid/llquid  interface.  The  analysis  is  possible  for  nonlinear,  as  well  as  linear  kinetics. 
It  is  found  that  linear,  parabolic,  and  exponential  kinetic  laws  should  yield  distinguishable 
waveforms  ot  interface  motion.  The  consistency  of  experimental  results  can  be  estab¬ 
lished  by  investigating  crystallization  kinetics  as  a  function  of  crystallographic  orienta¬ 
tion  and  of  square-wave  amplitude  and  fundamental  frequency. 


ACKNOWLEDGMENTS 

The  authors  are  indebted  to  Dr.  John  Simmons,  National  Bureau  of  Standards,  for 
his  useful  suggestion  which  led  to  the  derivation  of  Eq.  (21),  to  the  Research  Computation 
Center  staff  of  NRL  for  their  assistance  with  programming  the  CDC  380C  computer  for 
numerical  solutions,  and  to  the  Materials  Sciences  Division  of  the  Advanced  Research 
Projects  Agency  for  partial  financial  support  under  ARPA  Order  418. 


NRL  REPORT  6637 


13 


REFERENCES 

1.  Wilson,  H.A.,  Phil.  Mag.  50:23»  (1900) 

2.  Hillig,  W.B.,  and  Turnbull,  D.,  J.  Chem.  Phys.  24:914  (1956) 

3.  Volmer,  M.,  and  Marder,  M.,  Z.  Physik.  Chem.  (Leipzig)  A  154:97  (1931) 

4.  Kramer,  J.J.,  and  Tiller,  W.A.,  J.  Chem.  Phys.  37:841  (1962) 

5.  Kramer,  J.J.,  and  Tiller,  W.A.,  J.  Chem.  Phys.  42:257  (1965) 

6.  Rigney,  D.A.,  and  Blakely,  J.M.,  Acta  Met.  14:1375  (1966) 

7.  James,  D.W.,  pp.  767-773  in  Crystal  Growth,  ed.,  H.S.  Peiser,  New  York: Per gamon, 
1967 

8.  Tiller,  W.A.,  Acta  Met.  14:1383  (1966) 

9.  Rigney,  D.A.,  and  Blakely,  J.M.,  Acta  Met.  14:1384  (1966) 

10.  Utech,  H.P.,  and  Early,  J.G.,  Acta  Met.  15:1238  (1967) 

11.  Sekerka,  R.F.,  J.  Chem.  Phys.  46:2341  (1967) 

12.  Jakob,  M.,  Heat  Transfer,  Vol.  1,  New  York:Wiley,  1949 

13.  O'Hara,  S  ,  Tarshis,  L.A.,  and  Tiller,  W.A.,  J.  Chem.  Phys.  46:280C  (1967) 

14.  Burton,  W.K.,  Cabrera,  N.,  and  Frank,  F.C.,  Phil.  Trans.  Roy.  Soc.  A243:299  (1951) 

15.  Glicksman,  M.E.,  and  Void,  C.L.,  Acta  Met.  15:1409  (1967) 


DOCUMENT  CONTROL  DATA  •  R  &  D 

St'flinU'  «>/  filfv,  tauly  «»f  .il>vfr<i<#  diiii  .irifioUifinn  nmsf  hi'  eiijerrd  when  #ft«*  overall  repot!  is  C  humified) 


originating  activity  f  Corporate  .i.iiftor)  ca.  RLf’OHT  security  Classitic  a  hon 

Naval  Research  Laboratory  Unclassified 

Washington,  D.C.  20390 


I  IU.PORT  TITLE 


MEASUREMENT  OF  ATOMIC  KINETICS  OF  SOLIDIFICATION  USING 
PELTIER  HEATING  AND  COOLING 


4  DESCRIPTIVE  NOTES  ( Type  ot  report  and  inclusive  dares) 

This  is  a  final  report  on  one  phase  of  the  problem;  work  on  other  phases  is  continuing. 


S  AUTHORiSi  ff'H.il  niimr,  middle  initial,  fast  name) 

R.J.  Schaefer  and  M.E.  Glicksman 


December  5,  1967 


ea.  CONTRACT  OR  GRAN#  NO 


NRL  Problem  63M01-10 


ft.  PRQJEC  T  NO 


RR  007-01-46-5408 


ARPA  418 


70.  TOTAL  NO  OF  PAGES 


I  9a.  ORIGINATOR’S  REPORT  NUMHERIS) 


NRL  Report  6637 


I  9ft.  OTHER  REPORT  NOIS)  (Any  other  number*  that  may  be  assigned 
this  report) 


10  DISTRIBUTION  STATEMENT 


This  document  has  been  approved  for  public  release  and  sale; 
its  distribution  is  unlimited. 


(Office  of  Na1 
Research  Pr< 
in.  D.C.  2036 


13  abstract 


Periodic  Peltier  heating  and  cooling  may  be  generated  at.  a  solid/liquid  inter¬ 
face  by  the  passage  of  an  alternating  electric  current.  Melting  and  freezing  will 
then  occur  at  rates  which  wi’l  be  influenced  by  the  kinetic  law  of  molecular  attach¬ 
ment  relating  solidification  velocity  to  interface  temperature.  The  equation 
describing  the  interface  attachment  kinetics  may  be  determined  by  measuring  the 
response  of  a  solid/liquid  interface  to  square -wave  Peltier  heating  and  cooling. 

The  mathematical  methods  required  for  extracting  the  kinetic  equation  by  analysis 
of  the  waveform  of  interface  motior  v  that  this  technique  should  give  distinguish- 
ably  different  waveforms  for  linear,  t  i-abolic,  and  exponential  kinetics.  There  are 
advantages  of  the  method,  but  precautions  must  be  taken  to  eliminate  possible  con¬ 
fusing  effectr. 
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Solid/liquid  interface 
Interface  temperature 
Solidification  velocity 
Interface  attachment  kinetics 
Linear  freezing  kinetics 
Parabolic  freezing  kinetics 
Exponential  freezing  kinetics 


DD  ,"“.1473 


(BACK) 


(PAGE-  2) 


R  P  O  9  ,14-  2  39 


Security  Classifies!  ir 


